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The temperature dependence of the London penetration depth A was measured for an untwined 
single crystal of YBa2Cu307_^ along the three principal crystallographic directions (a, b, and c). 
Both in-plane components (X a and A&) show an inflection point in their temperature dependence 
which is absent in the component along the c direction (A c ). The data provide convincing evidence 
that the in-plane superconducting order parameter is a mixture of s + d— wave symmetry whereas 
it is exclusively s— wave along the c direction. In conjunction with previous results it is concluded 
that coupled s + d— order parameters are universal and intrinsic to cuprate superconductors. 

PACS numbers: 76.75.-K, 74.72.Dn, 74.25.Ha 



It is believed that the CuC>2 planes are the essen- 
tial building units in cuprate high-temperature supercon- 
ductors (HTS's) where superconductivity occurs. Even 
though either static or dynamic distortions of these 
planes destroy the cubic symmetry, many theoretical ap- 
proaches ignore the observed orthorhombicity and ideal- 
ize the planar structure, mostly in order to justify a pure 
d— wave order parameter. Early on it was, however, em- 
phasized that cuprates must have a more complex order 
parameter than just c?-wave [2, supported by many 
experiments like nuclear magnetic resonance (NMR) [3(, 
Raman scattering , angle-resolved electron tunneling 
Q, Andreev reflection 0], angular-resolved photoemis- 
sion (ARPES) 3, muon-spin rotation fuSR) HE El, 
and neutron crystal- field spectroscopy [12j . In addition, 



experiments along the c axis like, 



tunneling \U 



bi-crystal twist Josephson junctions T4(, optical pulsed 
probe , and optical reflectivity [lj| suggest that a pure 
s— wave order parameter is realized here. 

Multiple order parameter scenarios were proposed 
shortly after the BCS theory in order to account for a 
complex band structure and interband scattering (ItI 
181 Il9j |. This approach has the advantage that high- 
temperature superconductivity can easily be realized 
even within weak coupling theories since interband scat- 
tering provides a pairwise exchange between different 
bands which strongly enhances the transition tempera- 
ture (T c ) as compared to a single band model. The first 
realization of two-band superconductivity has been made 
in Nb doped SrTiC-3 and has not attracted very much 
attention. With the discovery of high-temperature super- 
conductivity in MgB2 two-gap superconductivity became 
more prominent, and meanwhile many more systems ex- 
hibiting multi-band superconductivity have been discov- 



ered (see, e.g., Refs. I20|,l2l|,l22|,l23|,l24|). Interestingly. 



in all these systems the coupled superconducting order 
parameters are of the same symmetry, i.e., s + s, d + d. 



In this respect HTS's are novel since here mixed order 
parameter symmetries are realized, namely, s + d. Theo- 
retically it has been shown that mixed order parameters 
support even high values of T c and lead to an almost 
doubling of the transition temperature as compared to 
coupled order parameters of the same symmetry [25| . 

In order to prove that complex order parameters are 
intrinsic and universal to HTS, previous /iSR measure- 
ments H, [l(], [H| were continued for another HTS fam- 
ily, namely YBa2Cu307_,5. The /iSR technique has 
the advantage that it is bulk sensitive and a direct 
probe of the London penetration depth which is highly 
anisotropic in HTS's. Recent results for Lai-zSr^CuO^ 
and YBa2Cu408 HE3, II | clearly demonstrate the exis- 
tence of two coupled s+d— wave gaps in the Cu02 planes 
and an s— wave gap along the c axis in YBa2Cu40s- 
While these findings already suggest that a complex 
gap structure is intrinsic to HTS, the new results on 
YBa2Cu307_<5, presented below, support this conclusion 
consistently. We are thus reasoning that s + d— wave 
superconductivity in the planes and s— wave supercon- 
ductivity along the c direction are intrinsic and universal 
to this complex material class. In addition, this finding 
imposes serious constraints for theoretical models, e.g., 
neither purely electronic nor 2D based approaches cap- 
ture this complicated picture. 

The crystal was grown by a crystal pulling technique 
[26| and exhibited a rectangular shape of approximate 
size of 4x4x1 mm 3 . The sample was detwinned by anneal- 
ing it under stress for 2 months at 400°C. The total frac- 
tion, where a and b axis are exchanged, occupies approx- 
imately 8 to 10% of the entire crystal, as confirmed by 
measurements with a polarized microscope. T c and the 
transition width were determined by DC-magnetization 
measurements and found to be 91.2 K and 2 K, respec- 
tively. 

The transverse-field //SR experiments were carried out 



2 



at the 7rM3 beam line at the Paul Schcrrer Institute 
(Villigcn, Switzerland). The samples were field cooled 
from above T c to 1.7 K in magnetic fields ranging from 
0.012 T to 0.64 T. The typical counting statistics were 
^20-25 million muon detections per experimental point. 
The experimental data were analyzed within the same 
scheme as described in Refs. This is based 

on a four component Gaussian fit of the /iSR time spec- 
tra where one component describes the background sig- 
nal stemming from muons stopped outside the sample, 
and the three other components describe the asymmet- 
ric local magnetic field P(B) distribution in the super- 
conductor in the mixed state. The magnetic field pen- 
etration depth A was derived from the second moment 
of P(B) since A~ 4 cx (AB 2 ) cx a 2 c The supercon- 

ducting part of the square root of the second moment 
(cr sc cx A -2 ) was obtained by subtracting the normal 
state nuclear moment contribution (onm) from the mea- 



sured a, as a 2 c = a 2 — a 2 m (see Ref. |9j for details) 
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FIG. 1: (Color online) Temperature dependences of a a 



b oc 



a; 



of YBa2Cu307-<5 measured after field cooling in 



HoH =0.1 T (a) and 0.64 T (b). The red lines represent 
results of a numerical calculation using the two-gap model 
S Ei EH with parameters as summarized in Table [I] The 
contributions of the small s— wave and the large d— wave gap 
to the in-plane superfluid density are shown by the blue and 
the black lines, respectively. 

Since cuprates are highly anisotropic, the relation 
a sc oc A~ 2 has to be extended to account for magnetic 



fields applied along the three crystallographic directions 
(i,j,k = a,b,c). For the field applied along the i— th 
principal axis the penetration depth is determined from 
the second moment like A~ fc 2 = (XjXk)^ 1 oc ajk = yjo] o% 
J28|. Here the index sc was omitted for clarity. The 
magnetic field dependence of the in-plane penetration 
depth has first been measured for different fields (0.05 T, 
0.1 T, 0.2 T, and 0.64 T). The field was applied along 
the crystallographic c axis and, subsequently, the sam- 
ple was cooled down from above T c to 1.7 K. For this 
field configuration the in-plane component of the super- 
fluid density a a b oc A~ b 2 = (AaAk)" 1 was obtained. In 
Fig. [1] a a b is shown as a function of temperature for 
two representative fields of 0.1 T and 0.64 T. In anal- 
ogy to previous results 0, E3, 11 1, an inflection point 
in cr a b(T) is observed at T ~ 10 K, which is a typical 
signature for the coexistence of a small s— wave and a 
large d— wave gap. Accordingly, the analysis of the data 
was performed by assuming that a{T) can be decom- 
posed into two components having d— wave and s— wave 

<? d {T) @, M 



11|. The 



symmetry as a(T) = er s (T) 
temperature dependences of the individual components 
were obtained within the same framework as presented in 
Refs.1, [ToL [ill . The comparison of experimental and the- 
oretical results is made in Fig.[TJ where the red lines refer 
to the sum of the two components, whereas the blue and 
the black lines display the individual s— and d— wave con- 
tributions, respectively. For all magnetic fields the anal- 
ysis was based on common zero-temperature gap values 
(Aq and Aq) but field dependent second moments [c s (0), 
a d (0)]. The zero-temperature gap values obtained in this 
way arc Ag = 0.707(11) mcV, A$ = 22.92(9) meV, and 
are in good agreement with results from tunneling ex- 
periments for the d— wave gap (see, e.g., Ref. [29h . The 
parameters obtained from the analysis are summarized 
in Table [J 

The d— wave contribution to the total superfluid den- 
sity co = <7 d (0)/[er s (0) + cr d (0)] increases with increasing 
field (see Fig. [2j) , similar to the field dependence observed 
for Lai.83Sro.i7Cu04 Q. This dependence can be under- 
stood by the fact that superconductivity is suppressed 
stronger in the s— wave band with increasing field than 
in the d— wave band Q. 

The individual components of the inverse squared pen- 
etration depth (A^ 2 , A^ 2 , and A^T 2 ) as functions of tem- 
perature were obtained by applying the magnetic field 
along the three crystallographic axes (0.012 T along a 
and b, and 0.1 T along c). This yields the axis-related 
superfluid densities according to [11( : 



Vik/vjk oc Aj 



(1) 



The results are shown in Fig. [3] (a). Obviously a a and 
Ob have a very similar temperature dependence, in par- 
ticular, the inflection point at T ~ 10 K and the lin- 
ear increase in an intermediate range of temperatures 
(60 K> T > 10 K). The temperature dependence of 
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FIG. 2: (Color online) The d— wave contribution to the in- 
plane superfluid density ui = a d (0)/[a B (0) + u d (0)] as a func- 



tion of the magnetic field of YBa2Cu307_«. 
guide to the eye. 
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TABLE I: Summary of the two-gap analysis for untwined 
single-crystal YBa2Cu307_i for the magnetic field applied 
along the c direction. The meaning of the parameters is - 
fioH: external magnetic field, u d (0) and u s (0): d— wave and 
s— wave contribution to the zero-temperature /iSR relaxation 
rate tr(0), u = a d (0)/[a s (0) + a d (0)]: the contribution of the 
large d— wave gap to the total in-plane superfluid density, A d : 
d— wave gap at T = K, Aq: s— wave gap at T = K. 



HoH a d {0) a s {0) 
(T) (MS" 1 ) (ms- 1 ) 



A d 
(meV) 



(meV) 



0.05 1.78(2) 4.80(7) 0.729(12) 

0.1 2.01(2) 5.53(6) 0.734(11) 

0.2 1.87(2) 5.63(7) 0.751(13) 22.92(9) 0.707(11) 

0.64 1.49(2) 5.01(7) 0.771(16) 



o~ c is qualitatively very different from the one of o~ a 
and <Th. Here a saturation is observed at T < 40 K. 
a a {T) and <Jb(T) can be well described by the two- 
component approach mentioned above where the same 
zero-temperature gap values were used. From this anal- 
ysis the following individual contributions from the s— 
and the d— wave components along the a and b axis are 
obtained: cr^(0) = 1.19 /zs _1 , <'(0) = 3.83 ^s" 1 and 
of (0) = 2.51 ,us~\ of (0) = 5 - 95 M 6 *" 1 - Since the relative 
contributions of the large d— wave component are almost 
the same along a and b directions, namely, u> a = 0.70, 
LOb = 0.76, it is plausible to conclude that not the CuO 
chains are the cause of the two-component behavior [3(| , 
but that this is an intrinsic property of cuprates. The 
same conclusions were reached from different experimen- 
tal techniques as, e.g., NMR Q, Raman scattering [1,0], 
and ARPES Q. In particular, Masui et al. Q showed 
that the s+d symmetry of the order parameter is required 



in order to describe the Raman data, even for tetragonal 
Tl 2 Ba 2 Cu0 6+( 5 HTS's. 
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3: (Color online) (a) Temperature dependences of a a oc 



(Tb oc Ai , and a c oc A c of YBa2Cu307 



obtained from 



cr(T) measured along the crystallographic a, b, and c direc- 
tions by using Eq. (fl]). Lines represent results of the analy- 
sis within the two-component (black and red lines) and one- 
component (blue line) models [lj|. (b) Temperature depen- 
dences of the anisot ropy p arameters y^, y ca , and 7 c t, obtained 
as -fij = Xi/Xj — \J (Tj/cTi (see text for details). 

The temperature dependence of the c axis related su- 
perfluid density a sc cx X~ 2 resembles very much the one 
observed for YBa 2 Cu40 8 [ll[ and follows closely the 
one expected for a single s— wave gap. The full blue 
curve in Fig. [3] (a) corresponds to an analysis based on 
an isotropic s— wave gap with Aq = 17.52 meV and 
cr c (0) = 0.183 /is -1 . The present results are in good 
agreement with previous findings of tunnelling experi- 
ments, where an s— wave gap along the c direction was 
reported [l3j |. The s— wave component along the c axis is 
not easily detectable by most experimental methods be- 
cause either very well oriented films should be prepared 
or bulk methods have to be used. Due to the fact that 
many experiments are surface sensitive only, and ab ori- 
ented samples and films are hardly available, these tech- 
niques are unable to see the s— wave component along 
the c axis. Its observation is, however, important since 
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the coupling of a major d— wave component in the ab 
plane to the s— wave component along the c axis mixes 
both symmetries in the planes already. In addition, the- 
ory predicts for this scenario, that a pure d— wave order 
parameter is never observable. On the other hand, also 
along the c direction an admixture of the d— wave order 
parameter has to be present [25|. This latter statement 
could provide an explanation for the optical conductivity 
spectra along the c direction where strongly anisotropic 
gap like features have been observed together with a fi- 
nite density of states at the Fermi energy 

Finally, the anisotropy along all three crystallographic 
directions is addressed. This can be calculated by using 
Eq. (TT]) and defining the anisotropy as 7.^ = = 
yjojjoi. The results are shown in Fig. [3] (b). While the 
in-plane anisotropy (7^ ) is almost constant for all tem- 
peratures, both out-of-plane components (pf ca and 7 c h) 
exhibit a sharp increase, as expected from Fig. [3] (a). In 
addition, it is seen that while 7 a t is almost always close to 
1.2, 7 ca , 7 C (, arc substantially larger and reach values be- 
tween 5 and 6.5 in the low-T limit. This high anisotropy 
is in very good agreement with previously reported values 
[33l . [34l ] . It reflects the fact that the CuC>2 planes have 
nearly Fermi liquid-like metallic properties whereas along 
the c direction mostly insulating behavior is observed. 

Our conclusions from the above presented data are 
manifold. Since s + d— wave symmetries of the super- 
conducting order parameter were observed previously in 
various cuprate families by various different techniques 
[I II !, U 0, S, M, EES El, the new ^SR data to- 
gether with earlier results on structurally different com- 
pounds 0, [l(J 11 1 support the idea that this behavior is 
intrinsic and universal. Similarly, the observation of an 
s— wave order parameter along the crystallographic c axis 
is proposed to be intrinsic as well. Specifically, this latter 
point emphasizes the importance of the third dimension 
for HTS's which was neglected in most of the theoretical 
models. In this context it is worth mentioning that the 
importance of the c axis has already been emphasized 
from ab-initio band structure calculations, where trends 
in T c were correlated with Cu02 apical oxygen distances 
[35| . Also, from first principles doping dependent com- 
putations of ARPES intensities, it was concluded that 
contributions from the c axis are of crucial importance in 
understanding the physics of HTS's [36]. On the other 
hand, the observation of mixed order parameters and 
more specifically, the additional s— wave component, re- 
quire that the lattice must be considered in the physics 
of HTS's. 
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